1 



Microscopic Reaction Dynamics at SPS and RHIC 

§ : Steffen A. Bass a 

a Department of Physics, Duke University, Durham, NC 27708-0305 
RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY11973 

< 

1. Transport Theory at RHIC 

Transport Theory offers the unique capability of connecting experimentally observable 
quantities in a relativistic heavy ion collision with its time evolution and reaction dynam- 
ics, thus giving crucial insights into the possible formation of a transient deconfined phase 
of hot and dense matter, the Quark-Gluon-Plasma (for reviews on QGP signatures, see 
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e-g- BE)- 

Figure [I] provides an overview of different transport theoretical ansatzes currently on 
the market for the description of a relativistic heavy ion collision at RHIC energies. The 
timeline shows a best case scenario for what to expect: the formation of a QGP with 
subsequent hadronization and freeze-out. Bands with solid lines denote the safe range of 
applicability for the respective transport approach, whereas dashed/dotted bands refer 
to areas in which the approach is still applied but where the assumptions on which the 
^ \ approach is based upon may be questionable or not valid anymore. 

The initial state and early pre- equilibrium phase are best described in Classical Yang- 
Mills theory (CYM) || or Lattice Gauge Transport (LGT) M calculations - only these 
classes of models treat the coherence of the initial state correctly, but do not provide any 
meaningful dynamics for the later reaction stages. 

The Parton Cascade Model (PCM) || and related pQCD approaches treat the initial 
state as incoherent parton configuration, but are very well suited for the pre- equilibrium 
phase and subsequent thermalization, leading to a QGP and hydrodynamic expansion. 
Microscopic degrees of freedom in this ansatz are quarks and gluons which are propagated 
according to a Boltzmann Equation with a collision term using leading order pQCD cross 
sections. Augmented with a cluster hadronization ansatz the PCM is applicable up to 
hadronization. The range of this model can be even further extended if it is combined 
with a hadronic cascade which treats the hadronic phase and freeze- out correctly J/J. 

Nuclear Fluid Dynamics (NFD, see e.g. ||||) is ideally suited for the QGP and hy- 
drodynamic expansion reaction phase, but breaks down in the later, dilute, stages of the 
reaction when the mean free paths of the hadrons become large and flavor degrees of 
freedom are important. The reach of NFD can also be extended by combining it with a 
microscopic hadronic cascade model - this kind of hybrid approach (dubbed hydro plus 
micro) was pioneered in |10| and has been now also taken up by other groups fLl| . It is 
to date the most successful approach for describing the soft physics at RHIC. The biggest 
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Figure 1. Transport theory approaches for RHIC and their range of applicability. 



advantage of NFD is that it directly incorporates an equation of state as input - one of 
it's largest disadvantages is that it requires thermalized initial conditions and one is not 
able to do an ab-initio calculation. 

Last but not least string and hadronic transport models [0,0 have been very suc- 
cessful in the AGS and SPS energy domains. They treat the early reaction phase as a 
superposition of hadron-hadron strings and are thus ill suited to describe the microscopic 
reaction dynamics of deconfined degrees of freedom. In the later reaction stages, however, 
they are the best suited approach, since they incorporate the full spectrum of degrees of 
freedom of a hadron gas (including flavor dependent cross sections). 

It is thus important to note that there is not a single transport theoretical ansatz 
currently available, which is able to cover the entire time-evolution of a collision at RHIC 
in one self-consistent approach. 

2. Kinetic Evolution 

The initial reaction dynamics are best described by partonic degrees of freedom. Fig- 
ure ^| shows the time-evolution of the average squared momentum transfer (Q 2 ) in parton- 
parton collisions for central Au+Au reactions at RHIC, based on a PCM calculation ]Tjj ]. 
Early rescatterings of incident partons occur at rather large values of Q 2 ~ 7 GeV 2 . 
However, between the time t c , m , ~ 0.5 fm/c and t c , m , ~ 4.5 fm/c the average (Q 2 ) de- 
creases rapidly to a value of about 4.5 GeV 2 until hadronization sets in. The strong 
time-dependence of (Q 2 ) shows that there is no single Q 2 scale which can be chosen to 
unambiguously set the pQCD scale in an ultra-relativistic heavy-ion collision. The re- 
action dynamics in the partonic phase is dominated by gluon-gluon scattering 43%) 
followed gluon fusion 27%) and quark-gluon scattering 26%). Quark-quark and 
quark- ant iquark scattering only contribute on a few % level. These numbers are fairly 
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Figure 2. Time evolution of parton scat- 
tering (Q 2 ) for central Au+Au collisions at 
RHIC in the Parton Cascade Model 



CD 



1.0 



0.5 



0.0 



K + /K Balance Function 




PYTHIA 

I I L 



A A i AAAAA 







1 



2 



Ay 



Figure 3. K + /K Balance Functions for 
T c = 165 MeV, T c = 225 MeV and in a 
hadronizing string picture 



independent of the system size and change only marginally between proton-proton and 
Au+Au reactions. 

While there is no doubt about the importance of partonic degrees of freedom for the 
initial, early, reaction stages, the crucial question is, how long the deconfmed state actu- 
ally existed and whether this time-span is long enough for thermalization and collective 
effects to occur. Balance functions offer a unique model-independent formalism to probe 
the time-scales of a deconfmed phase and subsequent hadronization |I5| . Late-stage pro- 
duction of quarks could be attributed to three mechanisms: formation of hadrons from 
gluons, conversion of the non-perturbative vacuum energy into particles, or hadronization 
of a quark gas at constant temperature. Hadronization of a quark gas should approxi- 
mately conserve the net number of particles due to the constraint of entropy conserva- 
tion. Since hadrons are formed of two or more quarks, creation of quark-antiquark pairs 
should accompany hadronization. All three mechanisms for late-stage quark production 
involve a change in the degrees of freedom. Therefore, any signal that pinpoints the time 
where quarks first appear in a collision would provide valuable insight into understanding 
whether a novel state of matter has been formed and persisted for a substantial time. 

The link between balance functions and the time at which quarks are created has a 
simple physical explanation. Charge- ant icharge pairs are created at the same location in 
space-time, and are correlated in rapidity due to the strong collective expansion inherent 
to a relativistic heavy ion collision. Pairs created earlier can separate further in rapidity 
due to the higher initial temperature and due to the diffusive interactions with other 
particles. The balance function, which describes the momentum of the accompanying 
antiparticle, quantifies this correlation. The balance function describes the conditional 
probability that a particle in the bin p± will be accompanied by a particle of opposite 
charge in the bin p^- 
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Figure 4. Multiplicity change in % due 
to hadronic rescattering for various hadron 
species at SPS and RHIC. The error-bars 
give an estimate of the systematic error. 
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Figure 5. Inverse hadron tjit slopes at 
Uc.m. — 0. The lines depict pure hydrody- 
namics whereas the symbols refer to data 
and hydro+micro calculations. 



B(p 2 \pi) = 2 M^kPi) - p{b,p 2 \b,p 1 ) + p(a,p 2 \b,p 1 ) - p(a,p 2 \a,p 1 )} , (1) 

where p(b,p 2 \a,p\) is the conditional probability of observing a particle of type b in bin 
p 2 given the existence of a particle of type a in bin p\. The label a might refer to all 
negative kaons with b referring to all positive kaons, or a might refer to all hadrons with 
a strange quark while b refers to all hadrons with an antistrange quark. 

Figure ||] shows K + /K~ balance functions as predicted in a simple Bjorken thermal 
model for two hadronization temperatures, 225 MeV and 165 MeV as well as for frag- 
menting strings (utilizing PYTHIA), which would represent the hadronization scenario 
of a hadronic/string picture similar to that of RQMD or UrQMD. Since particles from 
cooler systems have smaller thermal velocities, they are more strongly correlated in ra- 
pidity and result in narrower balance functions. A strong sensitivity to the hadronization 
temperature and time can be clearly observed. 



3. Flavor Evolution 

A central issue for the characterization of the deconfmed phase on the basis of the mea- 
surable final state hadrons is how much the hadronic composition of the system changes 
due to hadronic rescattering: figure |] shows the relative change (in %) of the multiplicity 
for various hadron species at SPS and RHIC between hadronization and freeze-out. As 
is to be expected, the state of rapid expansion prevailing at hadronization does not allow 
chemical equilibrium to hold down to much lower temperatures. The hadronic rescatter- 
ing changes the multiplicities by less than a factor of two, cf. also | l6f . Thus, we have 



first evidence that a QGP expanding and hadronizing as an ideal fluid produces a too 
rapidly expanding background for a hadron-fluid with known elementary cross- sect ions 
to maintain chemical equilibrium down to much lower temperatures than Tc- 
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Figure 6. Freeze-out transverse radius 
and time distribution d 2 N/ (rTdrT(ltf r ) for 
pions at RHIC. 
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Figure 7. .R out /-R S ide for RHIC initial condi- 
tions, function of at freeze-out (sym- 
bols) and at hadronization (lines). 



However, a closer look provides more insight into the chemical composition. The 
changes are most pronounced at the SPS, were the baryon-antibaryon asymmetry is high- 
est (since the net-baryon density at mid-rapidity is highest). This manifests e.g. in a 
reduction of the antiproton multiplicity by 40-50% due to baryon-antibaryon annihila- 
tion. A and 5 are affected in similar fashion. It remains a matter of debate whether these 
changes in the hadronic yields are consistent with the assumption of chemical freeze-out 
at the phase-boundary. 

However, despite the large amount of rescattering occuring in the hadronic phase, some 
hadron species may yet be sensitive to conditions at hadronization: figure ^] displays the 
inverse slope parameters T obtained by an exponential fit to dNi/cPrriTdy in the range 
rriT — Tfii < 1 GeV for SPS and RHIC in a hybrid hydro+micro calculation and compares 



them to SPS data ||17|| . The trend of the data, namely the "softer" spectra of S's and fi's 
as compared to a linear T(m) relation is reproduced reasonably well. This is in contrast 
to "pure" hydrodynamics with kinetic freeze-out on a common hypersurface (e.g. the 
T = 130 MeV isotherm), where the stiffness of the spectra increases linearly with mass as 
denoted by the lines in fig. [5| When going from SPS to RHIC energy, the model discussed 
here generally yields only a slight increase of the inverse slopes, although the specific 
entropy is larger by a factor of 4-5 ! The reason for this behavior is the first-order phase 
transition that softens the transverse expansion considerably. 

The softening of the spectra is caused by the hadron gas emerging from the hadroniza- 
tion of the QGP being almost "transparent" for the multiple strange baryons. This can 
be seen by calculating the average number of collisions different hadron species suffer in 
the hadronic phase: whereas JTs suffer on average only one hadronic interaction, iV's and 
A's suffer approx. 5-6 collisions with other hadrons before they freeze-out. 

Thus, one may conclude that the spectra of S's and especially Q's are practically unaf- 
fected by the hadronic reaction stage and closely resemble those on the phase boundary. 
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They therefore act as probes of the QGP expansion prior to hadronization and can be 
used to measure the expansion rate of the deconfined phase. 



4. Freeze-out 

Kinetic freeze-out, at which all (elastic) interactions cease, terminates the evolution 
of the reaction discussed in figure [l]. However, the freeze-out time of the system is ill- 
defined, since interactions cease locally on a particle per particle basis. Fig. ^| shows 
the distribution of freeze-out points of pions in the forward light-cone, as well as the 
hadronization hypersurface from where all hadrons emerge in a hydro+micro approach. 
Freeze-out is seen to occur in a four- dimensional region within the forward light-cone 
rather than on a three-dimensional "hypersurface" in space-time. Similar results have 
also been obtained within other microscopic transport models |18| when the initial state 
was not a QGP. It is clear that the hadronic system disintegrates slowly (as compared to 
e.g. the hadronization time), rather than emitting a "flash" of hadrons (predominantly 
pions) in an instantaneous decay. 

A first order phase transition leads to a prolonged hadronization time as compared to a 
cross-over or ideal hadron gas with no phase transition, and has been related to unusually 
large Hanbury-Brown-Twiss (HBT) radii |T9|-pl|]. The phase of coexisting hadrons and 
QGP reduces the "explosivity" of the high-density matter before hadronization, extending 
the emission duration of pions j|T9|-pT|]. This phenomenon should then depend on the 



hadronization (critical) temperature T c and the latent heat of the transition. For recent 
reviews on this topic we refer to p2p5f . 

It has been suggested that the ratio R on t/R S ide should increase strongly once the initial 
entropy density becomes substantially larger than that of the hadronic gas at T c [pill . 
The strong T c dependence of R ut/R S ide in such a purely hydrodynamical scenario can be 
seen in fig. |7| (solid and dotted lines). Here, the -R ou t/-Rside ratio at hadronization (calcu- 
lated in a hydrodynamical scenario) is compared to the ratio after subsequent hadronic 
rescattering and freeze-out in a combined hydro+micro approach |p4}| . Clearly, up to 
Kt ~ 200 MeV R OVLt / -Rside is independent of T c , if hadronic rescatterings are taken into 
account. Moreover, at higher Kt the dependence on T c is even reversed: for high T c the 
-R out /-R S ide ratio even exceeds that for low T c . Higher T c speeds up hadronization but on 
the other hand prolongs the dissipative hadronic phase that dominates the HBT radii. 
This is because during the non-ideal hadronic expansion the scale of spatial homogeneity 
of the pion density distribution increases, as the pions fly away from the center, but the 
transverse flow can hardly increase to counteract. Therefore, after hadronic rescattering 
i? out /-R S ide does not drop towards higher K T (in the range K T < 3m,r). 

For central collisions of Au nuclei at ^Js = 130A GeV, preliminary data from STAR 
gives Ront/Rsidc — 1-1 at small Kt P5|| . Results from RHIC at higher Kt will test 



whether a long-lived hadronic soft-rescattering stage, associated with the formation and 
hadronization of an equilibrated QGP state is indeed seen in heavy ion collisions at the 
highest presently attainable energies. 
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the results reviewed in this article. This work was supported by DOE grant DE-FG02- 
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